State Orthogonalization by Building a Hilbert Space: A New Approach to Electronic 

Quantum Transport in Molecular Wires* 
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Quantum descriptions of many complex systems are formulated most naturally in bases of states 
that are not mutually orthogonal. We introduce a general and powerful yet simple approach that 
facilitates solving such models exactly by embedding the non-orthogonal states in a new Hilbert space 
in which they are by definition mutually orthogonal. This novel approach is applied to electronic 
transport in molecular quantum wires and is used to predict conductance antiresonances of a new 
type that arise solely out of the non-orthogonality of the local orbitals on different sites of the wire. 
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The predictions of quantum mechanics that relate to 
observable phenomena do not depend on the particular 
basis that is selected to represent state vectors in Hilbert 
space. However choosing a basis of states whose phys- 
ical significance is clear and that are mutually orthogo- 
nal can be extremely helpful in theoretical work. Bases 
constructed from the eigenstates of a set of commuting 
operators that represent physical observables have both 
of these desirable properties [fij. But for complex sys- 
tems that are composed of simpler building blocks, it is 
tempting to use the eigenstates of the Hamiltonians of 
the separate building blocks as basis states even though 
such bases are not orthogonal. In solid state physics and 
quantum chemistry where the building blocks are atoms, 
each with its own electronic eigenstates, this choice of 
basis gives rise to the widely used tight-binding j|] and 
Hiickel || models. An analogous approach that mod- 
els nucleons as bags of quarks j| is used in theoretical 
work on nuclei and nuclear matter. |5| The representa- 
tions of the Hilbert spaces of complex systems that are 
obtained in this way are intuitively appealing but the 
non-orthogonality has been a significant drawback. Stan- 
dard orthogonalization schemes such as Gram-Schmidt 
do not help here because they are unwieldy for large sys- 
tems and do not preserve the atomistic character of the 
basis states. Wannier functions provide orthogonal- 
ized local basis states for perfect periodic structures, but 
they have the disadvantage of not being eigenstates of 
atomic Hamiltonians and also have no analog for disor- 
dered solids, liquids or molecules. Lowdin functions || do 
not require a periodic lattice but they too are not eigen- 
states of atomic Hamiltonians and do not have a simple 
physical interpretation. Thus, rather than working with 
nonorthogonal bases, pj it has been customary in much 
of the literature, for the sake of simplicity, to neglect the 
overlaps between the non-orthogonal states, as is done 



in LCAO (linear combination of atomic orbitals) mod- 
els of electronic structure, tight-binding theories of An- 
derson localization in disordered systems, and Hubbard 
and t-J models of electronic correlations. However, non- 
orthogonality can have non-trivial physical implications. 
For example, it is closely related to gauge interactions 
and fractional exclusion statistics, as has been pointed 
out by Haldane || . Better ways to treat it are therefore 
of interest. 

In this Letter we introduce a simple, powerful and gen- 
eral method that facilitates obtaining exact solutions of 
quantum problems in which the non-orthogonality of the 
basis is important. Our approach is to build a new Hilbert 
space around the non-orthogonal basis states, a space in 
which these states are by definition orthogonal, and to 
work in this new Hilbert space. We demonstrate the 
power, simplicity and flexibility of this novel approach 
by applying it to electronic quantum transport in molec- 
ular wires and predicting that these should exhibit con- 
ductance antiresonances of a new type. These antireso- 
nances are a direct and surprising physical consequence 
of the non-orthogonality of the electronic orbitals on dif- 
ferent atoms of the wire, which is included naturally in 
our theory. 

A molecular wire is a single molecule connecting a 
pair of metallic contacts. Such devices have recently 
begun to be realized experimentally and their electrical 
conductances are being measured. [pHl2[ Electron trans- 
port in molecular wires has been studied theoretically by 
considering the transmission probability for electrons to 
scatter through the structure. Jl0| , ^3| -p7t As with other 
mesoscopic systems, the electrical conductance G of the 
molecule is related to the transmission probability T at 
the Fermi level by the Landauer formula G = %-T. |l§| ] 
Molecular wires display a number of interesting trans- 
port phenomena, one of which is transmission antireso- 
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nances. ]l9|- pl| , [l4| Antiresonances are zeroes of the trans- 
mission T and correspond to electrons being perfectly 
reflected by the molecule. They also occur in semicon- 
ductor systems. 0] An analytic theory of molecular an- 
tiresonances was first proposed by Ratner Jl9[| in the con- 
text of electron transfer between donor and acceptor sites 
of a molecule, but his treatment did not include the ef- 
fects of the non-orthogonality of the atomic states. The 
present approach yields an analytic description of molec- 
ular wire transport (including antiresonances) that treats 
this non-orthogonality exactly and permits us to explore 
its physical implications. 

Our starting point is the Schrodinger equation 



H\V) = E\V) 



(1) 



where the Hamiltonian operator H and eigenvector \^S) 
are defined in a Hilbert space A. We wish to solve Eq. ([I]) 
for \^f) which is expressed as = ^ ^n\n) in a non- 
orthogonal basis {\n)} of Hilbert space A. In this basis 
Eq.(pl) takes the form 



(2) 



where H m „ = (m\H \n) and S m . n = (m\n) is the overlap 
matrix. |23[ Let us rewrite Eq.(0) as 



^H E ^ n = E% 



where 



n m,,n — n m.n 



m.n u m,n J 



(3) 



(4) 



Eq.(||) can be viewed as the matrix form of the 
Schrodinger equation 



H E \%') = E\V') 



(5) 



where H E and \^>') are a new Hamiltonian operator and 
its eigenvector denned in a new Hilbert space A' |24| ] 
in which the basis states {\n}} are orthonormal with 
(m\n) = 5 m . n , 1^') = En nas the same coefficients 

^ n as I and the new Hamiltonian operator is defined 
by its matrix elements through (m\H E \n) — H E n . Ac- 
cording to Eq.(|]), H E is Hermitian in A' because E is 
real and H m ^ n , S m . n and 5 m<n are Hermitian matrices. 

Thus we have transformed a problem that was formu- 
lated in terms of a nonorthogonal basis into an equivalent 
one in an orthogonal basis in a different Hilbert space. 
Other orthogonalization schemes (such as that of Lowdin 
H) differ from ours in this as well as other respects. j2J§] 

It should be noted that only the eigenvectors of 
that have the eigenvalue E have the same coefficients 
as eigenvectors of the true Hamiltonian H . The other 
eigenvectors of H E do not correspond to any eigenstate 
of the physical Hamiltonian H , but they never the less 
play an important role when calculating the Green's func- 
tion corresponding to H E . 



Since no assumptions at all have been made about the 
nature of the system being considered, this method of 
orthogonalization by switching to a new Hilbert space 
is extremely general. If the basis states {\n}} are tight 
binding atomic orbitals, then the present transformation 
(unlike the transformation to Wannier functions) can be 
used irrespective of the types of atoms involved or their 
locations in space. Furthermore, our transformation has 
the additional flexibility that the non-orthogonal basis 
states need not all be of the same generic type. For ex- 
ample, some of them may be atomic orbitals and others 
molecular orbitals on some cluster (s) of atoms that form 
a part of the physical system. This flexibility will be ex- 
ploited below. We now proceed to outline the application 
to molecular wire quantum transport. 

We begin by solving analytically an idealized model of 
a molecular wire consisting of a molecule attached to two 
identical semi-infinite single-channel leads which are rep- 
resented by ID chains of atoms. For this system we find 
it convenient to choose a non-orthogonal basis consisting 
of atomic orbitals {\n)} with n = -oo... - 1 on the left 
lead and n = 1 . . . oo on the right lead and molecular 
orbitals (MO's) {\j)} for the molecule. In terms of this 
basis we write electron eigenstates |\E r ) of H in which the 
electron is incident on the molecule from the left lead and 
transmitted with probability T through the molecule to 
the right lead as 



(6) 



The transmission probability T that enters the Landaucr 
electrical conductance of the wire is given by T — 'til 2 . 
Solving for analytically in the nonorthogonal basis 
is difficult so we transform to the new Hilbert space A' 
where the solution is more straightforward. As we have 
already shown, the coefficients remain the same in A' 
where the basis {\n), \ j)} is orthogonal, so that the above 
expression for T is valid in either Hilbert space. 

We evaluate "fi by solving a Lippmann-Schwinger (LS) 
equation that describes electron scattering in the molec- 
ular wire. We define this LS equation in the new Hilbert 
space A 1 where it takes the form 



|*') = |*'> +G'(E)W E \V). 



(7) 



Here is the eigenstate of the transformed Hamilto- 
nian H E whose coefficients correspond to those of the 
untransformed scattering eigenstate ys>) defined above 
Eq.(^|). W E is defined by separating H E into two parts, 
H E = H E + W E , where the matrix elements of H E be- 
tween lead orbitals {|n}} and MO's {\j)} all vanish in the 
space A' and W E couples the molecule to the adjacent 
lead sites. |<i>') is an eigenstate of H E with eigenvalue E 
that is confined to the left lead. G'(E) = (E - H^)' 1 is 
the Green's function of the decoupled system. 

The validity of the LS Eq. (|7j) depends crucially on the 
clear distinction between states on the leads and those on 
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the molecule that can only be made in Hilbert space A'; 
non-orthogonality leads to contradictions if analogs of the 
entities that enter (^) are constructed in A. The trans- 
formation to A' also introduces energy dependent hop- 
ping into the transformed Hamiltonian as prescribed in 
Eq. (|J). The energy dependence of the coupling W E be- 
tween the molecule and leads in Hilbert space A' will be 
important in the determination of antiresonances. Our 
choice of a set of MO's {\j)} that are mutually orthogo- 
nal in A means that the Green's function for the isolated 
molecule is formally unaffected by our transformation. 
This choice allows a simpler evaluation of the molecular 
Green's function that enters G'(E). An atomic orbital 
basis set non-orthogonal in A could be chosen instead, 
however the solution would then be less transparent. 

The evaluation of the Green's function G'(E) and so- 
lution of the LS equation ([?]) in the orthogonal basis 
{\ri), \j)} of Hilbert space A' is straightforward and will 
be presented in detail elsewhere. Here we will focus on 
the salient results for molecular wire conductance antires- 
onances and their physical significance. We find 

p & i 

*i = — (8) 

1 [(l-Q)(l-R)-PS] {> 

where 

3 
3 

R = G' 1A Y / (W\ J ) 2 G' J 

3 

S = G[ A J2w\ j G' j W^ 

3 

The sum over j is over only the MO's. In the above, 
Wf j = Hi j — ESij is the energy-dependent hopping el- 
ement of H E between the first lead site and the j th MO 
in terms of the hopping element of the original Hamil- 
tonian H and the overlap in the non-orthogonal basis. 
The Green's function on the molecule is expanded in 
terms of its molecular eigenstates (which, as mentioned 
above, are unchanged by the transformation) and this 
gives G'j = 1/(E — Cj) for the j th MO with energy ej. 
G[ 1 is the diagonal matrix element of the Green's func- 
tion G'(E) at the end site of the isolated lead. 

Conductance antiresonances of the molecular wire oc- 
cur where the transmission T = \^i\ 2 is equal to zero. 
From Eq. this happens when P = 0, i.e., at Fermi 
energies E that are the roots of 

y~* ~ ES i,j)( H j-i - ES 3-i) _ q ^ 
j ^ £j ' 

Two distinct mechanisms for antiresonances in molec- 
ular wire transport can be identified from Eq. (fj): 



In the first of these mechanisms, antiresonances arise 
due to an interference between molecular states that may 
differ in energy, as is seen directly from Eq. (^): An elec- 
tron incident from the left lead, hops from the lead site 
adjacent to the molecule onto each of the MO's with a 
weight W E _ 1 . It propagates through each of the dif- 
ferent orbitals j and hops onto the right lead with a 
weight Wfj . These processes interfere with each other 
and where they cancel ([)]) is satisfied and an antireso- 
nance occurs. This is in essence the same interference 
mechanism as has been identified previously in work on 
electron transfer between molecular donor and acceptor 
sites and, in the absence of the overlaps Sij and Sj-,—1, 
our result (||) agrees with that obtained there. jl9| 

The second antiresonance mechanism, which has no 
analog in previous work, arises solely from the nonorthog- 
onality of atomic orbitals that we have included in our 
analytic theory with the help of the Hilbert space trans- 
formation. It occurs when only a single MO k couples 
appreciably to the leads. In such cases, Eq. (^) becomes 
(H k ,-i - ESk,-i)(H hk - ESi, k ) = 0. Two antireso- 
nances are possible in this case. They occur at energies E 
where a matrix element of the transformed Hamiltonian 
H E — H — E(S — 1) that is responsible for hopping be- 
tween the molecule and one of the leads vanishes in the 
new Hilbert space. Thus the non-orthogonality of two 
orbitals can actually prevent electron hopping between 
these orbitals from taking place, blocking electron trans- 
mission along the wire and creating an antiresonance. 
This is a counterintuitive effect since one would normally 
expect orbital overlap to aid electron transfer between 
the orbitals rather than hinder it. It underscores the im- 
portance of including the effects of nonorthogonality fully 
in tight-binding theories. 

The above analytic theory of antiresonances was de- 
veloped for an idealized molecular wire model with semi- 
infinite single channel leads. We now compare these ana- 
lytic results with numerical calculations for a more real- 
istic molecular wire model. The system we consider con- 
sists of (100) Au leads bonded to a molecule as shown in 
the inset to Fig. . It is representative of a class of current 
experimental devices which use a mechanically controlled 
break junction to form a pair of nanoscale metallic con- 
tacts which are then bridged by a single molecule, the 
molecular wire. || The molecular wire we consider con- 
sists of two "chain" segments and an "active" segment. 
The purpose of the chains is to reduce the many propa- 
gating electron modes in the metallic contacts down to a 
single mode which propagates along the chains. Thus 
the (finite) chains supplant the ID ideal leads of our 
analytic model. We model this molecular wire and its 
bonding to the 3D metallic contacts by using extended 
Hiickel to calculate the hopping elements and overlaps 
between the non-orthogonal atomic orbitals that make 
up this system. It is the interaction between the non- 
orthogonal orbitals on the chains and the active segment 
that generates the antiresonances. Each chain consists 
of 7 C-H groups and is terminated with a sulphur atom 
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which bonds to a gold lead. For the energies of interest 
(near the Fermi level of gold) these chains only support 
a single ir mode. We chose an arbitrary active molec- 
ular segment with two 7r-like MO's which only interact 
with the 7r mode of the chains. The active segment is 
considered to be long enough that there is no direct cou- 
pling between the chains, as in our analytic theory. Fig. 
shows a plot of the contact-to-contact transmission cal- 
culated numerically for this model. The arrows indicate 
the locations of the antiresonances predicted by our an- 
alytic model (i.e., Eq. (^|)) using the same model param- 
eters. The agreement is very good; two antiresonances 
are found in each case at -10.2 eV and -10.4 eV, close 
to the Fermi energy of the gold leads. The transmission 
does not drop exactly to zero since in this calculation 
second nearest neighbor interactions are also included. 
The experimental signature is a drop in the differential 
conductance of the molecular wire. The agreement be- 
tween Eq. (||) and Fig. indicates that our analytic result 
derived for the idealized molecular wire model using the 
new approach to take account of non-orthogonality has 
predictive power for more complex systems. 

In conclusion: Many of the quantum problems that 
arise in physics and chemistry are formulated most nat- 
urally in a basis of states that are not mutually orthog- 
onal. In this Letter we have shown that the exact so- 
lution of such problems is greatly facilitated by embed- 
ding these non-orthogonal basis states in a new Hilbcrt 
space in which they are by definition mutually orthogonal 
but the matrix elements of the Hamiltonian are energy- 
dependent. The power, simplicity and flexibility of this 
novel approach was illustrated by applying it to analytic 
and numerical calculations of electronic quantum trans- 
port in molecular wires. A new mechanism for molecu- 
lar wire conductance antiresonances was identified which 
arises solely out of the non-orthogonality of local orbitals 
on different sites of the wire. 
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FIG. 1. Electronic transmission probability T calculated numerically for molecular wire shown in inset. Antiresonance is 
predicted at Al and A2 by the analytic condition (^|). The active molecule has two -k levels with energies e a — —13.0 eV and 
€b = —9.0 eV. The coupling and overlap of these two levels to the left (L) and right (R) chain molecules are Wh,a = —5.0 eV, 
W L ,b = -2.5 eV, S L ,a = 0.3, S L ,b = 0.2 and W R , a = -2.7 eV, W R , b = -1.8 eV, S R , a = 0.25, S R , b = 0.15. 
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